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ABSTRACT 

The  empirical  power  of  a  new  multivariate  goodness-of- 
fit  test  proposed  by  Foutz  (1980)  is  investigated.   The 
test  has  been  applied  to  Monte  Carlo  samples  from  bivariate 
and  trivariate  normal  distributions  with  a  variety  of  mean 
vectors  and  covariance  matrices.   The  null  hypothesis 
tested  is  that  the  sample  is  from  a  multivariate  normal 
distribution  with  0_  mean  vector  and  covariance  m.atrix  the 
identity  I_.   The  observed  number  of  rejections  in  5000 
replications  is  used  as  the  measure  of  effectiveness  of 
the  test.   The  results  indicate  that  the  Foutz  test  is 
quite  capable  of  detecting  mean  and  variance  shifts  but 
is  not  as  powerful  against  covariance  shifts. 
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I.   INTRODUCTION 

In  statistical  analysis,  choosing  the  correct  distri- 
bution to  model  available  data  is  of  importance.   A  class 
of  procedures  known  as  goodness-of-f it  tests  has  been 
derived  to  test  the  hypothesis  that  a  set  of  samples  is 
from  a  given  distribution.   Many  of  these  tests  are 
readily  available  and  are  well  known,  such  as  the  Chi- 
square  or  the  Kolmogorov-Smirnof f  (K.S.)  goodness-of-f it 
test.   These  tests  were  designed  for  univariate  distri- 
butions and  are  not  usable  as  multivariate  goodness-of-f it 
tests  in  their  present  form. 

In  1980  Robert  V.  Foutz  [Ref.  1]  proposed  a  new  multi- 
variate goodness-of-f it  test  that  will  be  called  the  Fn 
test  in  the  sequel.   In  analogy  to  the  K.S.  test  the  Fn 
test  compares  a  hypothesized  cumulative  distribution  func- 
tion (CDF)  with  a  "continuous  empirical  distribution 
function"  (CEDF)  formed  from  sampled  data.   Foutz  found 
the  null  distribution  of  the  test  to  be  distribution  free 
as  well  as  being  independent  of  the  number  of  variates  p. 

Foutz  obtained  an  integral  expression  for  the  null 
distribution  of  the  Fn  test  statistic  ,  and  closed  form 
solutions  for  sample  size  2  or  3  were  provided.   The 
complexity  of  the  integral  expression  increases  with 
sample  size,  and  a  normal  approximation  to  the  null  distri- 
bution was  given  for  use  with  larger  sample  sizes.   Although 
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the  Fn  test  was  designed  as  a  multivariate  goodness-of-f it 
test  it  can  also  be  used  to  fit  univariate  distributions. 
Franke  and  Jayachandran  [Ref.  2]  compared  the  empirical 
power  of  the  Fn  test  with  that  for  the  Chi-square  test  and 
the  K.S.  test.   The  results  indicated  that  the  Fn  test 
competes  well  with  these  other  tests. 

The  power  of  the  Fn  test  as  a  multivariate  goodness-of- 
fit  test  is  investigated  in  this  thesis.   A  description  of 
the  Foutz  test  is  given  in  Section  II  and  the  Monte  Carlo 
methods  of  simulation  are  presented  in  Section  III.   The 
results  and  conclusions  are  in  Section  IV.   A  Fortran  code 
for  the  application  of  the  Fn  test  is  available  in  the 
Appendix. 


II.   THE  FOUTZ  TEST 

The  Fn  test  for  multivariate  goodness-of-f it  is  based  on 
a  comparison  of  a  hypothesized  CDF  with  a  continuous  empiri- 
cal distribution  function  (CEDE)  derived  from  a  sample.   The 
first  step  in  the  determination  of  the  CEDE  is  the  construction 
of  what  are  known  as  statistically  equivalent  blocks.   A 
general  method  for  determining  statistically  equivalent  blocks, 
due  to  Anderson  [Ref.  3],  is  described  below. 

Given  a  random  sample  X.,Xr^,...,X   _-,  from  a  p-variate  con- 
tinuous distribution,  select  n-1  functions  h,  (X),  k  =  1,2,..., 
n-1,  not  necessarily  distinct,  such  that  each  h,  (X)  has  a 
continuous  distribution.   These  functions  are  referred  to  as 
cutting  functions  and  will  be  used  to  partition  the  sample 
space  into  blocks.   Let  k, ,k^,...,k   ,  be  a  permutation  of 

1,2,..., n-1.   Order  the  X. 's  according  to  h,  (X)  and  define 

1  K-L 

X(k,)  as  the  k, th  order  statistic.   The  sample  space  is  par- 
titioned into  two  blocks . 

B,   =     X:   h,   (X)   <   h,   (X(k,  )  ) 
±      I  -     ^1  -   -    ^1  -   ^   j 

B,   =    'X:   h,   (X)   >   h   (X(k,  )  )    . 
2       I  -     k^  -         T,^  _   ^    J 

At  the  second  step  if  0  <  k^  <  k   the  k-1  X's  in  B,  are 
ordered  according  to  h,  (X) ;  X(k2)  is  defined  as  the  k2th  in 
the  ordering.   Define  a  cut  on  B,  obtaining  3  blocks  as  follows: 
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^11   =   ^1  '     U--  \^^^    1^2^^^^^^^  ' 


i 


^12   =   ^  ^    ^''       \^^^     '    \^^^^2^^l    ' 


^20   ~   ^2* 


Now  consider  the  other  alternative,  k^  >  k-,  .   We  rank  the 
((n-l)-k-,)  X's  in  the  second  block  E^   according  to  h,  (X) 
and  let  X(k^)  be  the  (k--k, ) th  largest  in  the  ranking.   De- 
fining a  cut  at  h,  (X(k^))  we  obtain  the  3  blocks, 

^2  ~ 


^10   "   ^1' 


^21   =   ^2  ^    ^^  \    (X)  1^2  ^^^^2^^   ' 


y.   h  (Y\     >h   fyrvn' 


^22   =   ^2  '^  j=^=   \'^'  ^  \<^('^2"   • 

The  process  is  continued  until  all  the  cutting  functions  are 
exhausted.   This  results  in  a  partition  of  the  sample  space 
into  n  statistically  equivalent  blocks,  which  are  denoted  by 
B^,  i  =  1 , . . . ,n. 

In  the  univariate  case  an  intuitively  appealing  choice  for 
the  cutting  functions  is  the  identity  function  viz.,  h(X)  =  X 
for  all  k.   The  resulting  statistically  equivalent  blocks  are 
then  (-°°,X(1)]  ,  (X(l)  ,X(2)]  ,  .  .  .  ,  (X  (n-1)  ,+o°)   where  X(j)  is  the 
jth  order  statistic.   The  multivariate  analogue  is  to  choose 
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individual  coordinates  as  cutting  functions,  viz.,  h,  (X)  =  X  -^ 
the  jth  coordinate  of  X.   An  example  illustrating  the  con- 
struction of  the  blocks  in  the  bivariate  case  is  given  below 
for  a  sample  of  size  8. 

Let  (2,4,6,8,1,3,5,7)  be  the  permutation  vector  K.   Define 

h,  (X)  =  X    '  ^^^  first  coordinate  of  X,    for  k  =  2,4,6,8  and 

(2) 
h, (X)  =  X    ,  the  second  coordinate,  for  k  =  1,3,5,7.   Figure  1 

k  —    — 

gives  a  graphical  representation  of  the  rectangular  coordinate 
method  of  forming  blocks  and  Figure  2  is  the  representation 
for  the  polar  coordinate  method.   The  random  sample  that  was 
used  in  both  figures  is  found  in  Table  I. 


TABLE  I:   SAMPLE  BIVARIATE  DATA 

N  =  8  • 

Observation  12345678 

Coordinate 

1  -3.54       2.25      -1.00       .71      2.00    -    .75      -2.25       0.00 

2  0.00    -2.25         0.50       .00      1.25    -1.50      -1.50    -0.50 


The  first  element  of  the  permutation  vector  is  k  =  2  and 
h2 (X)  =  X   ,  therefore  xi    is  defined  to  be  the  second 

smallest  first  coordinate.   This  partitions  the  sample  space 
into  two  blocks. 


1       I  X.   X     1  Xj     , 

!  x=  x'^)  >  x'l' ! . 

(  j 


B    =   ,  X:   X'  '  >  X 
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FIGURE  1:   STATISTICALLY  EQUIVALENT  BLOCKS — 
RECTANGULAR  COORDINATES 
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FIGURE    2:       STATISTICALLY    EQUIVALENT    BLOCKS — 
POLAR   COORDINATES 
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The  second  element  of  the  permutation  vector  is  k^   =  4, 

h. (X)  =  X    and  k_  >  k, .   Hence  the  block  B-  is  partitioned 


into  two  sub-blocks. 


Sx:  X<"  <  X<"!  , 
(  ) 

jx:  X<^'  >  X<1)!  , 
(  ) 


^21   =   ^2  "  )^--       ^  -   ^2 


^22 


where  X^    is  the  second  largest  coordinate  among  the  X's  in 

block  B^.   At  this  stage  the  sample  space  is  partitioned  into 

three  blocks.   Next,  the  third  element  of  the  permutation 

vector  and  the  corresponding  cutting  function  define  another 

partition  of  one  of  the  three  blocks  into  two  sub-blocks. 

This  process  is  continued  until  the  permutation  vector  is 

exhausted,  at  which  stage  the  sample  space  will  be  partitioned 

into  9  statistically  equivalent  blocks. 

The  CEDF  is  now  constructed  by  spreading  a  mass  1/n  within 

each  block.   If  H„  is  the  hypothesized  CDF  and  H   the  CEDF, 

0  -^  "^  n 

the  test  statistic  Fn  takes  the  form 


Fn   =   SUP  lH„(X)  -  H„(X) | .  [1] 


Let  D.,  i  =  l,2,...,n,  be  the  probability  contents  of  the 

blocks  8.  under  the  null  hypothesis  H  ,  i.e.,  D.  =   /  dH^ ( x) 

i 

A  computational  form  of  the  Foutz  test  statistic  is. 
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n 


Fn   =    y   Max  (0,-  -  D. ) .  (2) 

i=l        n    1 


Foutz  gave  the  following  representation  for  the  cumulative 
distribution  of  the  test  statistic 

X  X 

P(Fn<x)   =    /   ...    /  g^{6^,S^,  .  .  .  ,6^_^)  dd-^dS^f  '  '  '  f^^^_i' 

— oo  — oo 

(3) 
where 


g„(6  w5^,  .  .  .  ,5^  n)   =   n!  (n-1)  : 
n   1   z       n— 1 


for 


k     >_     S^>    («2-^i)  >  ...  >  («n.i-^n.2)  >  "^n-l 


The  evaluation  of  this  integral  is  cumbersome  and  has  not  been 
carried  out  for  n  >  5 .   Foutz  has  therefore  derived  a  large 
sample  normal  approximation  given  by 


(1/2)  ,  -1, 

n   ^   (x  -  e   ) 

Lim  P[Fn<  x]   =   $[ ^^ TTTT^^  *         ^"^^ 

n-^°°  (2e    -  5e   ) 


To  check  the  accuracy  of  the  normal  approximation,  Franke 
and  Jayachandran  [Ref.  4]  generated  80,000  samples  of  sizes 
20,  30  and  50.   Table  II  contains  the  empirical  significance 
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TABLE  II:   EMPIRICAL  SIGNIFICANCE  LEVEL  OF  THE  FOUTZ  Fn  TEST 


Sample  Size 

Normal 
Significance 
Level 

.10 

.05 

.01 


20 


30 


50 


.0757 

.0800 

.0859 

.0372 

.0399 

.0428 

.0082 

.0083 

.0093 

levels,  when  the  normal  approximation  was  used  to  determine 
the  critical  values  for  the  Fn  test. 

It  is  clear  that  the  rejection  rates  given  in  Table  II 
are  consistently  lower  than  the  nominal  values.   More  accu- 
rate critical  values  were  therefore  determined  from  the 
80,000  Fn  values  and  are  presented  in  Table  III. 

TABLE  III:   APPROXIMATE  CRITICAL  VALUES  FOR  Fn  TEST 


Sample  Size 
Significance  Level 
.10 

.05 

.01 


20 


.42714 
(.43586) 

.44865 
(  .45513) 

.48659 
(.49127) 


30 


.41903 
(  .42383) 

.43553 
(  .43969) 

.46579 
(.46944) 


50 


.40816 
(.41150) 

.42116 
(  .42386) 

.44487 
(.44706) 


Values  in  parentheses  are  those  obtained  from  the  normal 
approximation  given  by  Foutz . 
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III.   DESCRIPTION  OF  THE  SIMULATION 

In  order  to  check  the  efficacy  of  the  Foutz  test  as  a 
multivariate  goodness-of-f it  test  a  simulation  was  run  to 
generate  sample  data  from  various  bivariate  and  trivariate 
normal  distributions.   The  hypothesis  tested  in  each  case 
is  that  the  sample  is  from  a  multivariate  normal  distribu- 
tion with  mean  vector  0^  and  covariance  matrix  the  identity 
I_.   Rectangular  and  the  polar/spherical  method  of  blocking 
were  both  used  and  compared  as  to  their  effect  in  each 
case. 

To  validate  the  blocking  schemes,  the  null  hypothesis 
is  tested  against  data  generated  from  the  distribution 
N(0_,I_).   Bivariate  and  trivariate  sample  sizes  of  20,  30 
and  50  are  used  to  compute  the  Fn  statistic  which  is  then 
compared  to  the  empirical  critical  levels  found  in  Table  III 
Rejection  rates  are  based  on  the  number  of  rejections  in 
20,000  replications  for  each  sample  size.   Comparing  the 
null  rejection  rates  to  the  nominal  significance  level 
used,  as  shown  in  Table  III,  provides  evidence  supporting 
both  blocking  methods  as  all  null  rejection  rates  are  close 
to  the  significance  level  used. 

The  empirical  power  of  the  test  was  then  investigated 
by  varying  the  distribution  tested.   This  investigation  is 
accomplished  in  three  different  ways.   First,  the  mean  is 
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shifted  away  from  the  £  vector  while  leaving  the  covariance  as 
the  identity  matrix.   This  is  done  to  investigate  the  ability 
of  the  test  to  detect  location  shifts.   The  covariance  matrix 
is  then  changed  from  the  identity  while  leaving  the  mean  as 
the  £  vector.   This  is  accomplished  by  changing  the  diagonal 
elements  alone  to  investigate  variance  shifts  and  then  shift- 
ing the  off  diagonal  elements  by  themselves  to  check  the  effect 
of  covariance  shifts.   A  primary  sample  size  of  20  was  chosen 
for  comparison  and  5000  replications  were  used  to  com.pute 
rejection  rates  for  each  distribution  tested.   Mixing  of  the 
three  types  of  shifts  is  also  simulated  to  investigate  the 
possible  confounding  effects  of  the  three  shifts.   Finally 
sample  sizes  of  30  and  50  are  run  on  a  few  of  the  distribu- 
tions to  determine  the  effect  of  increasing  the  sample  size. 

The  various  multivariate  normal  distributions  are  simulated 
in  the  following  manner.   Univariate  normal (0,1)  pseudorandom 
deviates  are  obtained  from  the  LLRZUSIDU  series  by  Lewis  [Ref .  5] 
and  grouped  to  form  a  multivariate  N  (£,!_)  p-variate  vector. 
Taking  the  }^  so  formed,  the  p-variate  N(0_,I_)  vector  random 
variable  is  transformed  by 


c"-^x*+  y  =  X  ,  (i; 


where 


C*  Z  C   =   I, 
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resulting  in  an  X  which  is  distributed  as  N(y_,E_)  .   The  Foutz 
test  is  then  applied  to  each  of  the  samples  consisting  of 
(n-1)  xs. 

An  example  using  a  bivariate  sample  helps  illustrate  the 
blocking  procedure  used.   Let  X^,Xy,...,X   _, ,  be  the  simu- 
lated bivariate  sample.   The  first  cut  is  made  on  X,    or 
the  first  coordinate  of  the  first  vector  X  .   Two  blocks  are 
formed, 


First  Second 

Coordinate      Coordinate 

(1) 


B-L   =   (-~,Xi  ^  (  —  ,+00) 


^2   "   ^-1^^  '^°°^  (-«',+<«)  . 


X_  is  taken  next  and  determined  to  be  contained  in  block  B, 

or  B- .   Suppose  X^  is  in  block  B-.   B-  is  then  partitioned 

(2) 
by  X-    or  the  second  coordinate  of  sample  X- .  Three  blocks 

are  now  defined  as, 


First  Second 

Coordinate      Coordinate 


^10   "   ^-°°'^{^^i         (  —  ,+00) 


^21   "   (^^^^+=°)         i-^'i^] 


B22   =   (2i{'^^'+"')         (2^2^,+-) 
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This  procedure  is  continued  by  examining  the  next  vector 
in  the  random  sample,  locating  the  block  that  it  is  contained 
in  and  partitioning  the  block  by  the  designated  coordinate. 
The  coordinate  cutting  functions  used  are  alternated  starting 
with  the  first  coordinate  for  the  first  cut.   Coordinate 
ranges,  as  shown,  are  used  to  designate  blocks  and  the  process 
is  continued  until  n  blocks  are  so  defined.   Given  any  random 
sample  this  method  can  be  shown  to  be  equivalent  to  a  unique 
permutation  vector  K  and  a  set  of  cutting  functions  (h,  }  as 
defined  in  Section  II. 

After  the  formation  of  the  statistically  equivalent  blocks, 
each  block  has  the  probability  content  of  1/n  and  must  be 
compared  to  the  hypothesized  content  using  the  statistic 


Fn   =    y   max  [0,--D.].  (2) 

.  ^,        n    1 
1=1 


D.,  the  probability  content  of  each  block,  under  the  null 
hypothesis,  is  defined  by  the  integral  of  the  null  density 
over  the  block.   The  integral  of  the  multivariate  normal  (0_,I_) 
over  a  rectangular  block  yields 


-P 

T^        -       f              f       f -,    \"^   -(l/2)x'Ix-, 
D.   =  j     ...    j        (2Tr)    e       dx. 

1 


(3) 
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This  reduces  to  the  product  of  the  marginal  densities  which 
may  be  easily  evaluated  with  many  available  routines,  elimi- 
nating the  need  for  numerical  integration. 

In  spherical  coordinates  D.  is  represented  by 


^2     ®  2     ^2  2 

D^   =    /      /      /   2tt^"^^^U^"^/^^P  sin((j))p^  dpd0d(|). 

H       H       pi 

(4) 


Upon  separation, 

^2  ^2     ,       ^2  ,  2  (-l/2)p^, 

D.   =     /    (1/2)  sin  4,  d4,    /  (271)"-^ d9    /    ^^-^ .   ,^      ^^ 

^      ({)]_  ej_  p^        (2tt)^^^ 

(5) 

Noting  that  with  a  change  of  variables  the  third  integrand 
is  a  Chi-square  density  with  3  degrees  of  freedom,  we  may 
use  a  closed  form  expression  to  evaluate  D  as  follows: 


D.   =   [|(cos  ,^2 -cos  4)i)]x[^(62-ei)]x[x3df(P2)-X3df(Pi)] 

(6) 


where 


X3df(Pi)   =   PtX3dflPi^'  i  =  1^2 


For  bivariate  data  the  use  of  polar  coordinates  leads  to 
similar  simplification  leaving  D.  in  the  form 
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°i   =  I?'92-6i)  .  (X2df(R2"  -  [X2df<Ri)l-        (7) 


After  the  calculation  of  the  probability  contents  D.  for 
the  n  blocks,  equation  (2)  is  used  to  evaluate  the  Fn  statis- 
tic for  each  generated  sample.   The  statistic  is  then  compared 
to  the  critical  values  found  in  Table  III  to  decide  if  the 
null  hypothesis  is  accepted  or  rejected.   Rejection  rates 
are  defined  by  the  number  of  rejections  divided  by  the  number 
of  replications  in  a  given  run.   The  rejection  rates  thereby 
define  an  empirical  power  for  the  simulated  distribution. 

The  major  component  of  the  Fortran  simulation  program 
used  to  evaluate  the  Foutz  statistic  for  a  given  sample  is 
available  in  the  Appendix.   It  has  been  adapted  for  use  for 
sample  sizes  up  to  50,  with  redimensioning  being  needed  for 
larger  sample  sizes.   The  progran  is  applicable  for  fitting 
data  from  any  hypothesized  multivariate  normal  distribution 
and  provides  the  Fn  statistic  as  computed  by  both  blocking 
methods  presented.   The  code  is  self-contained  except  for 
three  IMSL  routines,  LUDECP ,  MDNOR,  and  MDCH  [Ref .  6] .   These 
subroutines  provide  matrix  decomposition,  univariate  normal 
probabilities  and  chi-square  probabilities,  respectively, 
and  must  be  available  or  substituted  prior  to  utilization 
of  the  program. 
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IV.   RESULTS  AND  CONCLUSIONS 

The  results  of  the  simulation  are  summarized  in 
Tables  IV-XIV.   Rejection  rates  are  given  by  the  distribu- 
tion tested  and  the  significance  level  used.   Empirical 
power  curves  are  presented  in  Figures  3-8.   Rejection  rates 
are  plotted  against  the  magnitude  of  the  shift  in  mean, 
variance  and  covariance  for  the  distribution  tested.   All 
power  curves  are  based  on  5000  replicated  samples  and 
were  compared  at  the  a  =  .05  significance  level. 

The  results  for  the  case  in  which  the  distribution  of 
the  samples  is  the  same  as  the  hypothesized  distribution 
viz.,  N(0_,I_)  are  given  in  Tables  IV  and  V.   The  rejection 
levels  obtained  are  close  to  the  nominal  significance 
level  for  both  blocking  methods.   No  distinct  pattern  of 
variation  about  the  prescribed  levels  is  discernible  for 
either  method,  as  expected. 

The  rejection  rates  for  mean  shifts  are  given  in  Tables 
VI-VII  and  Figures  3-4.   Shifts  in  the  mean  vector  are 
detected  well;  a  shift  of  one  standard  deviation  in  a 
single  coordinate  resulted  in  a  60%  rejection  rate  for 
bivariate  or  trivariate  data.   Greater  shifts  in  mean  led 
to  even  higher  rejection  rates.   The  rectangular  method 
of  blocking  consistently  gave  about  a  10%  improvement  over 
the  polar/spherical  method  in  detecting  mean  shifts. 
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Results  for  variance  shifts  are  contained  in  Tables 
VIII  and  IX  and  the  power  curves  are  given  in  Figures  5 
and  6.   The  Foutz  test  did  not  detect  small  variance 
shifts  very  well  but  the  performance  of  the  test  was  far 
better  for  larger  shifts  or  shifts  in  more  than  one  coor- 
dinate.  No  one  method  of  blocking  performed  better  in  all 
cases  but  in  general  the  polar/spherical  method  seemed  to 
outperform  the  rectangular  method  for  detecting  variance 
shifts. 

The  results  for  changes  in  covariance  are  summarized 
in  Tables  X,  XI  and  Figure  7.   Covariance  shifts  are  not 
detected  well  for  either  blocking  method  except  for  highly 
correlated  data  with  the  correlation  coefficient  equal  to 
.9.   The  polar/spherical  coordinate  blocking  method  appeared 
to  perform  a  little  better  than  the  rectangular  coordinate 
method  of  blocking,  but  in  general  the  simulation  revealed 
that  the  Fn  test  is  not  very  powerful  against  covariance 
shifts . 

The  empirical  power  for  combinations  of  shifts  in  mean 
and  variance  or  covariance  are  presented  in  Tables  XII  and 
XIII.   Entries  are  based  on  an  a  =  .05  significance  level 
and  are  tabled  by  the  mean  vector  and  covariance  matrix 
of  the  sample  data.   Entries  farther  down  and  to  the  right 
correspond  to  greater  shifts  in  mean  and  variance/covariance 
and  are  generally  larger,  as  is  to  be  expected.   There  are 
no  apparent  confounding  problems  due  to  shifts  in  both 
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parameters.   The  rectangular  method  of  blocking,  however, 
did  outperform  the  polar/spherical  method  for  most  cases 
of  multiple  shifts. 

The  results  indicative  of  the  effect  of  increasing  the 
sample  size  are  summarized  in  Tables  XIV  and  XV.   Results 
for  sample  sizes  of  20,  30,  and  50  are  given  for  some 
representative  cases.   The  tables  reveal  higher  rejection 
rates  for  larger  sample  sizes  with  increases  being  compa- 
rable for  both  blocking  methods. 

This  study  was  limited  to  the  two  and  three  variate 
normal  distribution.   There  are  many  problems  for  further 
research.   Of  primary  concern  is  the  generation  of  percen- 
tage points  of  Fn  for  various  values  of  n.   The  intracta- 
■bility  of  the  problem  of  obtaining  the  exact  distribution 
requires  an  empirical  approach  to  finding  a  correction  to 
the  asymptotic  approximation  given  by  Foutz.   Since  the 
use  of  coordinates  as  cutting  functions  worked  well,  the 
method  should  be  tried  for  other  distributions  and  higher 
dimensions . 

In  conclusion,  the  Fn  test   is  found  to  be  a  viable 
option  for  testing  goodness-of-f it  of  multivariate  normal 
distributions.   These  encouraging  empirical  results  indicate 
further  study  should  be  conducted  to  explore  the  potential 
of  this  test  for  other  distributions. 
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TABLE  IV:   NULL  EMPIRICAL  REJECTION  LEVELS  FOR 
THE  BIVARIATE  NORMAL  DISTRIBUTION 


Significance  Level       .01  .05  .10 

Blocking  Method 

N  =  20 

Rectangular  .0098  .0488  .0940 

Polar  .0096  .0482  .1020 

N  =  30 

Rectangular  .0110  .0510  .0944 

Polar  .0082  .0454  .0890 

N  =  50 

Rectangular  .0120  .0498  .0950 

Polar  .0098  .0484  .0958 


BASED  ON  20,000  REPLICATIONS 
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TABLE  V:   NULL  EMPIRICAL  REJECTION  LEVELS  FOR 
THE  TRIVARIATE  NORMAL  DISTRIBUTION 


Significance  Level       .01  .05  .10 

Blocking  Method 


N  =  20 

Rectangular  .0104  .0440  .0982 

Spherical  .0120  .0518  .1048 

N  =  30 

Rectangular  .0114  .0480  .0956 

Spherical  .0140  .0484  .0914 

N  =  50 

Rectangular  .0098  .0484  .0960 

Spherical  .0088  .0478  .0914 


BASED  ON  2  0,000  REPLICATIONS 
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TABLE  VI:   REJECTION  RATES  FOR  SHIFTS  IN  MEAN  (BIVARIATE) 


N  =  20 


Critical  Value 

.01 

.05 

.10 

Mean  Tested 

0 

.0094 

.0488 

.0998 

0 

.0102 

.0482 

.0996 

-.5 

.0566 

.1684 

.2816 

0 

.0346 

.1096 

.2138 

.5 

.0574 

.1710 

.2700 

0 

.0430 

.1388 

.2298 

-.5 

.1408 

.3164 

.4534 

-.5 

.1038 

.2592 

.3610 

.5 

.1294 

.3024 

.4406 

.5 

.1230 

.3164 

.4534 

-1 

.4357 

.6664 

.7834 

0 

.2340 

.4484 

.6046 

1 

.4464 

.6700 

.7842 

0 

.2444 

.6664 

.7834 

-1 

.8382 

.9418 

.9748 

-1 

.7780 

.9212 

.9610 

1 

.8428 

.9418 

.9718 

1 

.6930 

.9212 

.9610 

-2 

.9936 

.9980 

.9996 

0 

.9926 

.9990 

.9996 

2 

.9948 

.9998 

1 

.0000 

0 

.9762 

.9950 

.9980 

-2 

1.0000 

1.0000 

1 

.0000 

-2 

1.0000 

1.0000 

1 

.0000 

2 

1.0000 

1.0000 

1 

.0000 

2 

1.0000 

1.0000 

1 

.0000 

BASED  ON  5000  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — POLAR 
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TABLE  VII;   REJECTION  RATES  FOR  SHIFTS  IN  iMEAN  (TRIVARIATE! 


N  =  20 


Critical  Value 

.01 

.05 

Mean  Tested 

0 

.0104 

.0440 

0 

.0120 

.0518 

0 

-.5 

.0492 

.1502 

0 

.0216 

.0980 

0 

.5 

.0480 

.1438 

0 

.0280 

.1036 

0 

-.5 

.1076 

.2704 

- .  5 

.0472 

.1516 

0 

.5 

.0972 

.2424 

.5 

.0658 

.1912 

0 

-.5 

.1826 

.3788 

-.5 

.0848 

.2198 

-.5 

.5 

.1738 

.3642 

.5 

.0848 

.2198 

.5 

-1 

.3782 

.5984 

0 

.1184 

.2942 

0 

1 

.3728 

.5984 

0 

.1174 

.2866 

0 

-1 

.7392 

.8892 

-1 

.3808 

.6020 

0 

1 

.7400 

.8892 

1 

.4670 

.6918 

0 

.10 


0982 
1048 

2474 
1632 

2484 
1874 

3990 

2516 

3688 
3046 

5170 
3584 

4948 
3584 

7212 
4212 

7212 
4234 

9410 
7338 

9410 
7934 
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TABLE  VII  (Continued) 


Critical 

Value 

.01 

.05 

.10 

Mean  Tested 

-2 

0 
0 

.9636 
.7772 

.9872 
.9138 

.9958 
.9916 

2 
0 
0 

.8992 
.8486 

.9676 
.9448 

.9832 
.9736 

-1 
-1 
-1 

.9134 
.7688 

.9778 
.8744 

.9894 
.9312 

1 
1 

1 

.9102 
.7936 

.9746 
.9244 

.9900 
.9598 

-2 
-2 

0 

1, 

.0000 
.9984 

1.0000 
.9996 

1, 
1, 

.0000 
.0000 

2 
2 

0 

1 

.0000 
.9998 

1.0000 
1.0000 

1, 
1, 

.0000 
.0000 

-2 
-2 
-2 

1 

1 

.0000 
.0000 

1.0000 
1.0000 

1, 
1, 

.0000 
.0000 

2 
2 
2 

1, 
1 

.0000 
.0000 

l.OOQO 
1.0000 

1. 
1, 

.0000 
.0000 

BASED  ON  5000  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — SPHERICAL 
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TABLE  VIII.   REJECTION  RATES  FOR  SHIFTS  IN 
VARIANCE  (BIVARIATE) 


N  =  20 


Critical  Values        .01  .05  .10 

Variance  Tested 

.0998 
.0996 

.5150 
.4640 

.5628 
.4928 

.7474 
.7574 

.8580 
.8368 

.8556 
.8856 

.9222 
.9600 

.9832 
.9804 

.9998 
.9988 


1 

0 

0 
1 

1 

0 

0 
3 

2 
0 

0 
2 

1 
0 

0 

5 

3 
0 

0 
3 

1 
0 

0 
7 

1 
0 

0 
10 

5 
0 

0 
5 

0 
0 

0 
10 

0094 

.0488 

0102 

.0482 

1864 

.3786 

1578 

.3342 

2228 

.4292 

1714 

.3582 

4030 

.6322 

4286 

.6448 

5790 

.7666 

5338 

.7450 

5640 

.7608 

6312 

.8106 

7092 

.8618 

8088 

.9228 

8998 

.9664 

8998 

.9665 

9956 

.9994 

9920 

.9978 

BASED  ON  50  0  0  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — POLAR 
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TABLE  IX.   REJECTION  RATES  FOR  SHIFTS  IN  VARIANCE 
(TRIVARIATE) 


N  =  20 


01  .05  .10 


0982 
1048 

3626 
4736 

4644 
3308 

5372 
7550 

6332 
8792 

6326 

5814 

7514 
6566 

9038 
8966 

7186 

9506 

9346 
9182 

9872 
9852 


Critical 

Value 

Variance 

Tested 

1 

0 

0 

0 

1 

0 

0 

0 

1 

3 

0 

0 

0 

1 

0 

0 

0 

1 

2 

0 

0 

0 

2 

0 

0 

0 

1 

5 

0 

0 

0 

1 

0 

0 

0 

1 

7 

0 

0 

0 

1 

0 

0 

0 

1 

2 

0 

0 

0 

2 

0 

0 

0 

2 

3 

0 

0 

0 

3 

0 

0 

0 

1 

3 

0 

0 

0 

3 

0 

0 

0 

3 

10 

0 

0 

0 

1 

0 

0 

0 

1 

5 

0 

0 

0 

5 

0 

0 

0 

1 

5 

0 

0 

0 

5 

0 

0 

0 

5 

0104 

.0440 

0120 

.0518 

0924 

.2372 

1606 

.3438 

1500 

.3330 

0888 

.2208 

1940 

.3832 

4146 

.6374 

2708 

.5026 

6100 

.8032 

2758 

.5012 

2510 

.4538 

4140 

.6270 

3240 

.5394 

6622 

.8372 

6752 

.8312 

3716 

.5980 

7880 

.9078 

7558 

.8896 

7256 

.8660 

9390 

.9770 

9292 

.9762 
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TABLE  IX   (Continued) 


N  =  20 


Critical 

Value 

.01 

Variance 

Tested 

10    0 
0   10 
0    0 

0 
0 
0 

.9572 
.9470 

10    0 
0   10 
0    0 

0 

0 

10 

.9972 
.9858 

.05  .10 


9866 

.9950 

9832 

.9926 

9998 

1.0000 

9970 

.9994 

BASED  ON  5000  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — SPHERICAL 
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TABLE  X:   REJECTION  RATES  FOR  SHIFTS  IN  COVARIANCE 
(BIVARIATE) 


N  =  20 


Critical  Value 

.01 

Covariance 

Tested 

1 
0 

0 
1 

.0094 
.0102 

1 
-.3 

-.3 

1 

.0152 
.0126 

1 
.3 

.3 

1 

.0126 
.0136 

1 
-.6 

-.6 
1 

.0288 
.0514 

1 
.6 

.6 

1 

.0250 
.0648 

1 
-.9 

-.9 

1 

.1166 
.2378 

1 
.9 

.9 

1 

.1122 
.2378 

.05  .10 


0488  .0998 

0482  .0996 

0558  .1068 

0598  .1274 

0576  .1178 

0656  .1258 

1008  .1782 

1576  .2560 

0912  .1702 

1838  .2984 

2996  .4446 

2982  .6162 

2996  .4446 

4710  .6042 


BASED  ON  5000  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — POLAR 
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TABLE  XI:   REJECTION  RATES  FOR  SHIFTS  IN  COVARIANCE 
(TRIVARIATE) 


N  =  20 

Critical 

Value 

.01 

.05 

.10 

Covariance  Tested 

1 
0 
0 

0 
1 
0 

0 
0 

1 

.0104 
.0120 

.0440 
.0518 

.0982 
.1048 

1 

0 

-.3 

0 
1 

-.3 

0 
1 

.0104 
.0124 

.0540 
.0488 

.1076 
.1066 

1 

0 

.3 

0 
1 
0 

.3 
0 
1 

.0106 
.0124 

.0468 
.0512 

.0972 
.1086 

1 
.3 
.3 

.3 

1 

.3 

.3 
.3 

1 

.0126 
.0162 

.0560 
.0676 

.1112 
.1292 

1 

0 

-.6 

0 
1 
0 

-.6 
0 

1 

.0152 
.0122 

.0740 
.0584 

.1394 
.1158 

1 

0 

.6 

0 

1 

0 

.6 
0 

1 

.0148 
.0128 

.0674 
.0582 

.1298 
.1194 

1 
.5 
.6 

.6 

1 

.6 

.6 

.6 

1 

.0308 
.0432 

.1136 
.1434 

.1960 
.2486 

1 
0 

-.9 

0 
1 

0 

-.9 

0 

1 

.0412 
.0254 

.1358 

.0974 

.2314 

.1842 

1 

0 

.9 

0 

1 

0 

.9 
0 

1 

.0402 
.0258 

.1386 
.1386 

.2368 
.2368 

1 
.9 
.9 

.9 

1 

.9 

.9 
.9 

1 

.1406 
.3646 

.3454 
.5950 

.4942 
.7122 

BASED  ON  50  00  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — SPHERICAL 
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TABLE  XII:   REJECTION  RATES  FOR  MULTIPLE  SHIFTS 
IN  MEAN  AND  VARIANCE-COVARIANCE 
(BIVARIATE) 


N  =  20 


Sigma 

1   0 
0   1 

1   .6 
.6    1 

2    0 
0    2 

2    .849 
.849     2 

5    1.34 
1.34    5 

Mean 

0 
0 

.0488 
.0482 

.0912 
.1176 

.1986 
.1522 

.2500 
.2162 

.9658 
.9572 

.5 
0 

.1710 
.1388 

.2398 
.2482 

.3110 
.2498 

.3702 
.3402 

.9720 
.9650 

1 
0 

.5606 
.4348 

.7346 
.5952 

.6384 
.5334 

.6828 
.6316 

.9820 
.9764 

1 
1 

.9418 
.8576 

.8774 
.8588 

.9350 
.8722 

.8658 
.3308 

.9892 
.9840 

2 
0 

.9998 
.9950 

.9998 
.9990 

.9902 
.9772 

.9950 
.9882 

.9990 
.9964 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — POLAR 

BASED  ON  5000  REPLICATIONS 
a  =  .05 
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TABLE  XIII: 


REJECTION  RATES  FOR  MULTIPLE  SHIFTS 
IN  MEAN  AND  VARIANCE -COVARIANCE 
(TRIVARIATE) 


Sigma 


10  0 
0  10 
0     0     1 


1     0     .6 

0     10 

,6     0       1 


N  =  20 

5  0  0 
0  10 
0     0     1 


10     0   .95  5     0     0 

0     10  0     5     0 

,95     0     1  0     0     5 


Mean 


0 
0 
0 

.0440 
.0518 

.0674 
.0582 

.5392 
.4584 

.7828 
.7840 

.9770 
.9720 

.5 
0 
0 

.0480 
.0280 

.1830 
.1176 

.5708 
.5034 

.7946 
.8020 

.9832 
.9740 

1 
0 
0 

.3728 
.1174 

.6352 
.2912 

.6852 

.6254 

.8206 
.8422 

.9888 
.9824 

1 
1 

0 

.7400 
.7392 

.9074 
.7454 

.9270 
.8602 

.9668 

.9454 

.9930 
.9870 

2 

0 

1 

.9982 
.9726 

.9956 
.9752 

.9716 
.9742 

.9736 
.9774 

.9978 
.9976 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY  —  SPHERICAL 

BASED  ON  5000  REPLICATIONS 
a  =  .  05 
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Sample 

size 

Shift 

.5 

0 

.5 

.5 

1   .3 

.3    1 

1    0 

0    3 

TABLE  XIV:   REJECTION  RATES  FOR  INCREASING 
SAMPLE  SIZES  (BIVARIATE) 


20  30  50 


a   =    .01 


0574 

.0860 

0430 

.0564 

1294 

.2026 

1230 

.1418 

0126 

.0140 

0136 

.0152 

1864 

.2722 

1578 

.2244 

2700 

.3228 

2298 

.2658 

4406 

.5424 

4534 

.4336 

1178 

.1174 

1258 

.1196 

5150 

.6132 

4640 

.5566 

1270 
0754 

3652 
2508 

0176 
0170 

4522 

3744 


a  =  .05 

.5  .1710  .2170  .2914 

0  .1388  .1630  .2238 

.5  .3024  .4144  .6030 

.5  .3164  .3076  .4826 

1   .3  .0576  .0624  .0728 

.3    1  .0656  .0624  .0760 

1    0  .3786  .4884  .6756 

0    3  .3342  .4304  .6016 

a  =  .10 

.5  .2700  .3228  .4256 

0  .2298  .2658  .3400 

.5  .4406  .5424  .7190 

.5  .4534  .4336  .6066 


1  .3  .1178  .1174  .1396 

.3  1  .1258  .1196  .1422 

1  0  .5150  .6132  .7800 

0  3  .4640  .5566  .7160 


BASED  ON  5000  REPLICATIONS 

FIRST  ENTRY — RECTANGULAR 
SECOND  ENTRY — POLAR 


44 


TABLE  XV:   REJECTION  RATES  FOR  INCREASING 
SAMPLE  SIZES  (TRIVARIATE) 


Sample  size 

Shift 

.5 

0 

0 

.5 

.5 

.5 

1 

0  .3 

0 

1   0 

.3 

0   1 

3 

0   0 

0 

1   0 

0 

0   1 

20 

30 

a  =  .01 

0480 

.0680 

0280 

.0362 

1738 

.2932 

0848 

.1662 

0106 

.0138 

0124 

.0134 

1606 

.2054 

0838 

.1138 

50 


1036 
0526 

5040 
3428 

0148 
0144 

3528 
2080 


.5 

0 

0 

.5 

.  5 

.5 

1 

0  . 

,3 

0 

1 

0 

3 

0 

1 

3 

0 

0 

0 

1 

0 

0 

0 

1 

a  =  .05 

1438  .1974  .2742 

1036  .1256  .1646 

3642  .5118  .7268 

2198  .3588  .5868 

0468  .0588  .0656 

0512  .0488  .0540 

3438  .3976  .5734 

2074  .2708  .4126 


.5 

0 

0 

.5 

.5 

.  5 

1 

0  . 

.3 

0 

1 

0 

.3 

0 

1 

3 

0 

0 

0 

1 

0 

0 

0 

1 

a  =  .10 

2484  .3024  .3876 

1874  .2132  .2650 


.4948  .6396  .8272 

.3584  .4912  .7040 

.0972  .1142  .1232 

.1086  .1030  .1102 

.4736  .5264  .6880 

.3156  .3880  .5450 


BASED  ON  50  00  REPLICATIONS   FIRST  ENTRY:  RECTANGULAR 

SECOND  ENTRY:   SPHERICAL 
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APPENDIX  A 
USER  REQUIREMENTS  AND  INPUT  FORMAT  FOR  PROGRAM  FOUTZ 

The  use  of  the  Computer  program  contained  in  Appendix  B 
requires  the  sample  size,  number  of  variates,  applicable 
data  and  the  Multivariate  Normal  distribution  being  tested 
as  described  by  the  mean  vector  and  the  variance-covariance 
matrix.   The  variables  containing  the  required  inputs  as 
well  as  the  required  input  format  are  as  shown  below. 

DESCRIPTION  OF  VARIABLES 

N Sample  size 

M Number  of  Variables  (2  or  3) 

SIGMAl Variance-covariance  Matrix 

(MxM) 

Bl Mean  Vector  (Mxl) 

X Matrix  of  Sample  Data  (MxN) 

INPUT  FORMAT 

N,M (215) 

SIGMAl (3F12.6)  Input  M  Rows 

Bl (F12.6)  Input  M  Rows 

X (3F12.6)  Input  M  Rows 

Input  data  is  echoed  in  the  output  providing  a  check 
for  correct  entry  of  data  as  well  as  is  the  decomposition 
of  SIGMAl.   The  Fn  statistic  as  computed  by  both  methods 
of  blocking  follows  completing  the  output  given  for  a 
single  run.   An  example  run  is  given  for  Trivariate  data 
of  sample  size  10. 
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SAMPLE  TRIVARIATE  RUN 


FOUTZ  TEST  FOR   3  VARIATE  NORJ'IAL 
THE  NUMBER  OF  OBSERVATIONS  =10 
OBSERVATIONS  ENTERED  AS  FOLLOWS 


3. 

.170000 

7, 

.540000 

4. 

.230000 

4. 

.160000 

5. 

.500000 

5. 

,580000 

2, 

.330000 

2. 

.910000 

6. 

.620000 

2. 

.530000 

3. 

.440000 

5. 

.660000 

1. 

.990000 

2. 

.630000 

6. 

,320000 

2. 

.260000 

2. 

.800000 

6. 

,730000 

2. 

.630000 

0. 

.290000 

6. 

.550000 

3. 

,440000 

4. 

.860000 

3. 

,150000 

3. 

,500000 

4. 

.670000 

8. 

,310000 

3. 

.580000 

3. 

.230000 

4. 

,970000 

R] 

[BUTION 

TESTED 

COVARIANCE  MATRIX 

1.000000  1.000000 
1.000000  3.000000 
1.000000     1.000000 

MEAN  VECTOR 
2.000000 
3.000000 
4.000000 


DECOMPOSITION 
1.000000 
-0.707107 
-0.500000 


OF  SIGMA 
0.0 

0.707107 
0.0 


1.000000 
1.000000 
5.000000 


0 
0 
500000 


WITH  POLAR  OR  SPHERICAL  COORDINATES 
FOUTZ  STAT=      0.593289 

WITH  RECTANGULAR  COORDINATES 
FOUTZ  STAT=      0.556877 
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c 


APPENDIX  B 
CCMPUTER  PROGRAM  FOUTZ 


C  ♦        FORTRAN  CODE  ♦ 

C  ♦  FCUT2  ♦ 

C 
C 
C 

C  N  SAMPLE    SIZE 

C  M  DIMENSION    OF    EACH    VECTOR 

C  IRAQ  N    VECTOR    DESIGNATING    COORDINATE    TO   CUT    ON 

C  SIGMAl  (M,M)COVARIANCE    MATRIX       TEST    DISTRIBUTION 

C  Bl  (M.IJMEAN    VECTOR 

C 

C 

C  MAIN    PROGRAM 

C 

C  PURPOSE  ! 

C  'reads     in    N,M    and    DIMENSIONS 

C  LIMITED    TO    N=50,M=3    AS    SET 


DIMENSION    IRAD{52)  ,VECT(50f6),WKVEC(6),8L0CK(51,12) 


$SIGMAl(3j3) ,Bl(3,i),X(50,3) ,TRAN(3,1),XTT(3, 1),C(3,3J, 
$BLCC{51,i2),XTTR(3a) 

REA0(5,990)N,M 
990  F0RMAT(2I5) 

NN=N+1 

MM=2*M 

NM1=N-1 

DO    10    1=1, N,M 

DO    5    J=1,M 

IRAD( J+I-1)=J 
5      CONTINUE 
10     CONTINUE 

CALL  DDRIVE(IRAD,VECT,WKVEC,8L0CK,BLCC,SIGMA1,B1,N,M, 
$NN»MMtTPAN,XTT,C,XTTR,X) 

STOP 

END 

C 

C  SUBROUTINE  DCRIVE 

C 

C  PURPOSEl 

C  DRIVES  PROGRAM  AND  VARIABLE  OIMENSINS  BASED  ON 

C  M  AND  N.   READS  IN  Bl ,  SIGMAl  AND  DATA  TO  BE 

C  TESTED.   ECHCS  INPUT  DATA  AND  PRINTS  THE 

C  RESULTING  FN  STATISTIC. 

C 

c 


SUBROUTINE    CDP  IVE(  IR AD , VECT , WKVEC, 8LCCK, oLOC , SIGMAl, 
$81 tN,M,NN,MMTTRAN,XTT,C,XTTR,Xi 

DIMENSION  IRACCN) , VECT ( N ,M ) , WKVEC( 3 J  , BLOCKC NN , 6) , 
$SIGMA1{M,M) ,Bi(M,l),TRAN(M,l ),X(N,M) , BLOC (NN , MM ) , 
$XTT(M, 1)  ,C(M,M),XTTR<M,  1) 

00    30    1=1, M 

REA0(5,992) (SIGMAl (I,J)  ,J=1,M) 
30  CONTINUE 

992         FGRMAT(3F12.6  ) 

DO    ^0    1=1, M 
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RE/iD(5, 993)81  (  1,1  ) 
993         F0RMAT(F12.6) 
40  CONTINUE 

DO    70    I=liN 

R6AD(5,995}(X(I, J)  ,J  =  1,M) 
70  CONTINUE 

CC  ECHO    INPUT    DATA 

WRITE(6,800)M 

WRITE(6»801)N 

800  FORMAT(  •!•  ,»FCLT2    TEST    FCR    'tlZ,*    VARIATE    NORMAL') 

801  FORMAT( 'O* ,'THE    NUMBER    CF    OBSERVATIONS    =«,I3) 

804  FORMATC     •,3F12.6) 

805  FORMAT! 'O*  , 'MEAN    VECTOR') 

806  FORMAT{»     'jFia.d) 
URITE(6t79t) 

791  FORMAT(  'OS 'OeSERVATIONS    ENTERED    AS    FCLLQWS»  ) 
CO    792     1  =  1  ,N 

792  WRITE(6,804){ X(I ,J) , J=1,M) 
WRITE(6,807) 

WRITE (6,808) 

807  FORMAT(  'O*  ,»DISTRIBUTICN    TESTED*) 

808  FORMAT(  'O', 'CCVARIANCE    MATRIX') 
DO    793    1=1, M 

793  WRITE<6Ta04){ SIGMA  1(1, J)  ,J=1,M) 
WRITE (6,805) 

WRITE  (  6,806)  (  (Bid  ,J)  ,J=  1,1)  ,1=1  ,M) 
CALL    OECOMP(SIGMA1,M,C,0) 
C  CALL    TRANSFCR^ATICN       ROUTINES 

DO    750    1=1, N 
DO    751    J=1,M 

751  XTT(J,1)=X( I, J) 

752  FORMAT(  •     •  ,3X,F12.6) 

CALL    TRANS(M,XTT,B1,TRAN,C,XTTR) 

DO    760    J=1,M 

VECT( I, J)=TRAN(J, 1) 
760         CONTINUE 
750         CONTINUE 

CALL    BLCCKS(VECT,N,NN,M,MM,IRAD,BLOCi<,WKVEC) 

CALL    3LCCKR(VECT,N,NN,M,MM,  IRAD,  BLOCWKVEC) 
C 
C    BLOCK    BY    POLAR    CR    SPHERICAL    COORDINATES 

CALL    FCUTZ(8LCCK,NN,MM,FN,MI 

WRITE(6,989) 

989  FORMAT( 'O* ,'WITH    PCLAR    CR    SPHERICAL    COORDINATES') 
WRITE(6,990)FN 

990  FORMATC     "j'FCUTZ    STAT=    •,F12.6) 
C    BLOCK    BY    RECTANGULAR    COORDINATES 

CALL    FCUTR(BLCC,NN,MM,FN,M) 
WRITE(6,988) 
988  FORMATC '0' , 'WITH    RECTANGULAR    COORDINATES') 

WRITE(6,99i)FN 

991  FORMATC     •,'FCUTZ    STAT=    ',F12.6) 
RETURN 

END 

C 

C  SUBROUTINE  CECOMP 

C 

C  PUR  PCS  E  I 

C  DECOMPOSES  THE  CCVARIANCE  MATRIX  ENTERED. 

C  USES  ChCLESKY  OECOMOSITION  VIA  IMSL  ROUTINE 

C  LUDECP  TO  PROVIDE  A  MATRIX  C  NEEDED  BY  THE 

C  ROUTINE  TRANS. 

C 

c 

SUBROUTINE    DECCHP ( SI GMA ,M, C, I V ) 

DIMENSION    SIGMA(M,M)  ,C ( M  ,M) • A( 51 ) ,UL (51 ) , LI ( 6 ) ,M 1 ( 6 ) 

IJ  =  1 
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DO    100    1=1, M 

CO    110    J=1,I 

A(  IJ)=SIGMA(I ,J) 

IJ=IJ+1 
110         CONTINUE 
100         CONTINUE 

CALL    LUDECP(A,UL,M,Dl,D2,ieR) 

00    120     1=1, M 

DO    130    J=1tM 

II=I*(I-1)/2>J 

IF(J.LT.I)C(I,J)=UL(II} 

IF<J.EQ,I)C(I  ,J)=1./UL{II) 

IF(J.GT.I)C( I,JI=0.0 
130         CONTINUE 

120  CONTINUE 
IFdV.EG.DGO    TO    121 
CALL     INVT(C,M,0,L1,M1) 
WRITE{6,799) 

799         FORMAT( 'O* ,»DECOMPCSITICN    OF    SIGMA" I 

DO    765    1=1, M 
765         WRITE(6,700)(C(I,J ), J=1,MJ 
700         FOPMATC     •,3F12.6) 

121  CONTINUE 
RETURN 
END 

C 

C  

C 

C  SUBROUTINE     INVT 

C 

C  PURPOSE 

C  INVERT    A    MATRIX 

C 

C  

c 


SUBROUTINE    INVT( A , N , D, L ,MI 
DIMENSION    A( 1),L(1),M( 1) 


c 

c 

SEARCH  FOR 

LARGEST  ELEMENT 

0=1.0 

NK=-N 

DO  80  K=1,N 

NK=NK+N 

L(K)=K 

M( K) =K 

KK=NK+K 

eiGA=A(KK) 

DO  20  J=K,N 

IZ=N«( J-l) 

DO  20  I=K,N 

IJ=IZ+I 

10 

IF(  ABS(BIGA)- 

ABS(A(IJ)))  15 

15 

BIGA=A(IJ) 
L(K)  =  I 
,M<K)=J 

c 
c 
c 

20 

CONTINUE 

INTERCHANGE 

ROWS 

J=L{K) 

IF(J-K)  35,35, 

25 

25 

KI=K-N 

CO  30  1=1, N 

KI=KH-N 

HOLD=-A{KI ) 

JI=KI-K+J 

20,20 
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A(KI )=A(JI) 
30  A(JI)  =HCLD 
C 

C         INTERCHANGE  COLUMNS 
C 

35  I=M(K) 

IF(I-K)     45,45,38 
38    JP=N*(I-1) 
DO    40    J=1,N 
JK=NK+J 
JI=JP■^J 
hOLD=-A{JK) 
A( JK)=A( JI) 
40    A(JI)    =HCLD 
C 

C  DIVIDE    COLUMN    BY    MINUS    PIVOT     (VALUE   OF    PIVOT 

C  ELEMENT    IS    CONTAINED    IN    BIGA) 

C 

45  IF(BIGA)    48,46,48 

46  C=0,0 
RETURN 

48    00    55    1=1, N 

IF(I-K)  50,55,50 
50    IK=NK+I 

A(  IK)=A(  lK)/(-BIGA) 
55    CONTINUE 
C 

C  REDUCE    MATRIX 

C 

DO    65    1=1, N 
IK=NK-H 
HOLD  =  A(  IK) 
IJ=I-N 
DO    65    J=1,N 
IJ=IJ+N 

IF(I-K)  60,65,60 
60  IF(J-K)  62,65,62 
62    KJ=IJ-I+K 

A( IJ)=HCLD*A(KJ)>A(IJ) 
65    CONTINUE 
C 

C  DIVIDE    ROW    EY    PIVCT 

C 

KJ=K-N 
DO    75    J=1,N 
KJ=KJ+N 

IF(J-K)  70,75,70 
70  A(KJ)=A(KJ)/BIGA 
75  CONTINUE 
C 

C         PRODUCT  OF  PIVCTS 
C 

D=D«BIGA 
C 

C         REPLACE  PIVOT  BY  RECIPROCAL 
C 

A(KK)=1.0/aiGA 
80  CONTINUE 
C 
C         FINAL  ROW  AND  COLUMN  INTERCHANGE 


C 


K=N 
100  K=(K-1) 

IF(K)  150,150,105 
105  I=L(K) 

IF(I-K)  120,120,108 
108  jg=N*(K-l) 
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JR=N*(I-1) 

00    110    J=1,N 

JK=JQ+J 

HaLO=A( JK) 

JI=JR+J 

^(JK)=-A(JI ) 
110    A(JI)     =HCLD 
120    J=M(K) 

IF(J-K)     100,100fl25 
125    KI=K-N 

00    130    1=1, N 

KI=KI+N 

HOLD=A(Kn 

JI=KI-K+J 

A(KI )=-A( JI) 
130    A(JI)     =HCLO 

GO    TO    100 
150    RETURN 

END 

C 

C  SUBROUTINE  TRANS 

C 

C      PURPOSE:  TO  TRANSFORM  OBSERVATIONS  TC  N(0,I) 

C  UNDER  THE  NULL  HYPOTHESIS.   USES  INPUT 

C  VALUES  OF  Bl  ANC  THE  MATRIX  C  FROM  DECOMP 

C  TO  TRANSFORM  THE  DATA  ENTERED  USING, 

C  * 

C  X   =  C(X-Bl). 

C 

c 


SUBROUTINE  TR ANS( M,XTT ,B1,TRAN,C ,XTTP ) 
DIMENSION  81(^^,1)  ,XTT{M,  1)  ,XTTR(M,1J  ,TRAN(M,1}  ,C(M,M) 
CALL  SUB{XTT,ei,XTTR,M,l} 
CALL  PRD(C,XTTR,TRAN,M,M,1) 
RETURN 
END 
C      

c 

C  SUBROUTINE  SUB 

C 

C  PURPCSE 

C  SUBTRACT  ONE  MATRIX  FROM  ANOTHER  TC 

C  FCRM  RESULTANT  MATRIX. 

C  

C 


SUBROUTINE    SU E { A , B ,R ,N,M ) 

DIMENSION    A(l)  ,8(1),R( 1) 
C 

C  CALCULATE    ^UMB£R    CF    ELEMENTS 

C 

NM=N*M 
C 

C         SUBTRACT  MATRICES 
C 

00  10  1=1, NM 
10  R(  I)  =  A(  I)-8(I) 

RETURN 

END 
C 
c  

c 

C  SUBROUTINE  PRO 

C 

C  PURPCSE 

C  MULTIPLY  TWO  MATRICES  TO  FORM  A 

C  RESULTANT  MATRIX. 

C 
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SUBROUTINE    PRC  (A  ♦  B  ,R  ,N  .M  ,L  ) 
DIMENSION    A(l)  ,B(l)fR( 1) 

IR=0 
IK=-M 

00    10    K=1,L 
IK=IK+M 
DO    10    J=1,N 
IR=IR+1 
JI=J-N 
I3  =  IK 
R(  IR)=0 
DC    10    1=1, M 
JI=JI+N 
18=13+1 
10    R(IR)=R(IR)+A(JI)*B(  IB) 
RETURN 
END 


C  SL8R0LTINE    BLOCKS 

C 

C  PURPOSES 

C  THIS  SUBROUTINE  TAKES  N  M-VARIATE  VECTORS  AND  PARTITIONS 

C  A  SPACE  OF  DIMENSION  M  INTO  N+1  STATISTICALLY  EQUIVALENT 

C  BLCCKS  BY  RECCRCING  BLOCK  COORDINATE  RANGES  IN  A  MATRIX 

C  BLOCK  BY  THE  USE  OF  SPHERICAL  OR  POLAR  COOROI NATESAS 

C  AS  CUTTING  FUNCTIONS.   THE  CUTTING  COORDINATE  USED  AT 

C  EACH  STEP  IS  CONTAINED  IN  A  VECTOR  IRAD. 

C 

c 

SUBROUTINE    BLCCKS ( VECT , N  ,NN ,M,MM  ,  IRA C  , BLOCK, WKVEC ) 
DIMENSION    VECT(N,M),6L0CK(NN,6),IRA0(N)  ,t^KVEC(6) 
ZL=1.0E-8 
BLCCKC 1, 1)=0.0 
eLOCK(l,2)=1000. 
BLOCKd  ,3J=G.O 
BLCCK( 1,4)=6. 2831853 
BLOCKd  ,5)=0.0 
BLOCKd  ,6)  =3.  1415927 
DO    100    J=1,N 
TEMP=0. 0 
CO    110    1=1, M 
110         TEMP  =  TEVP+VECT(J,n**2 
RAD=TEMP**.5 
IF  (RAD.GT.ZDGC    TO    112 
T0EG=6. 2831852 
P0EG  =  3.  1A15S27 
GO    TO    111 

112  TARG=VECT( J,l  )/( ( VECT( J , 1 ) **2  +  VECT( J , 2)**2)**  .  5) 
PARG=VECT(J,M)/RAD 

IF(TARG.LT.-1  .)TARG=-1 . 

IF{TARG.LT.l. ) GO    TO    1122 

TARG=1.0 

IF{PARG.LT.-1  .  )PARG=-1. 

1122  IF(PARG.LT.1.0)GO    TO    1123 
PARG=1.0 

1123  OEG=ACOS(TARG) 

IF(VECT( J,2).GT.0. JGO    TO    111 
IF(0EG.LT.1.5707963)G0    TO    113 
D£G=0EG+1. 5707963 
GO    TO    111 

113  DEG=0EG+4. 712389 

111   continue 

poeg=accs(parg; 

WKVECd  )=TEMP 
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WKVEC(2)=0EG 

UKVEC(3)=PDEG 

DO    120    1=1, NN 

IF(WKVEC(1  )  .GT.BLOCKdf  2))G0    TO    120 

IF(WKVEC(1).LT.BL0CK(I ,1))G0    TO    120 

IF(WKVEC(2) .LT.BLCCK{I,3) JGO    TO     120 

IF  (WKVEC(2J .GT.BLCCK{I,4)iG0    TO    120 

IF(i^.EQ.2)G0    TC    119 

IF  (WKVECO)  .LT.BLaCK(I,5)>G0    TO    120 

IF(WKVEC(3) .GT.BLCCK{I»6i)G0    TO    120 

119  CONTINUE 
IBLOCK^I 
GO    TO    150 

120  CONTINUE 
150  JJ=IRAD(J) 

BLIM=^*KVEC(  JJ) 

00    160     1=1, MM 
160         BLCCK( J+l,I)=eLQCK( IBLCCKfl ) 

BLOCK ( I  BLOCK, 2*JJi=BLIM 

BLCCK( J+1,2*JJ-1)=BLIM 
100         CONTINUE 

RETURN 

END 

C 

C  SUBROUTINE    FCUTZ 

C 

C  PURPOSE:    TO    CCMPUTE    THE    FOUTZ    STATISTIC    FROM    THE 

C  BLOCKS    DETERMINED    BY    SUBROUTINE    BLOCKS 

C  METHOD    USES    IMSL    ROUTINE    MOCh    TO    EVALUATE 

C  CHI-SGUARE    PROBABILITIES    TO    EVALUATE    THE 

C  CLOSED    FORM    EXPRESSION    GIVEN    FOR    0.    THE    FN 

C  STATISTIC     IS    GENERATED    BY    FCUTZ'S    CLOSED 

C  COMPUTATIONAL    FORMULA. 

C 

C 

SUBROUT  INE    FOUTZ ( BLOCK ,NN, MM, FN, M) 

DIMENSION    3L0CK(NN,6),P(51) 

CF=FLOAT(M) 

TP=0.0 

DO    100    1  =  1, NN 

CALL  MDCh(aLOCK( I,  1) ,0F,P1,IER) 

CALL  MDCH(BLQCK( 1,2)  ,DF,P2, lER) 

P3=P2-P1 

P4= (BLOCK (1, 4) -BLOCK ( 1 , 3 ) )/6 .283 1853 

IF(M.EQ.2JG0  TC  85 

P5= (COS ( BLOCK (  1,5) }-COS( BLOCK ( 1,6)) ) /2.0 

GO  TO  86 

85  P5=1.0 

86  P( I)=P3*P4*P5 
TP=TP+P( n 

100    CONTINUE 

FN=0.0 

DO  300  1=1, NN 

AMAX=1.0/NN-P(I) 

IF  ( AMAX.LT.O.)GO  TO  300 

FN=FN+AMAX 
300  CONTINUE 

RETURN 

END 

C 

C  SUBROUTINE  BLOCKR 

C 

C    PU  RPCSE I 

C  THIS    SUBROUTINE    TAKES    N    M-VARIATE    VECTORS    AND    PARTITIONS 

C  A    SPACE    OF    DIMENSION    M    INTO    N+1    STATISTICALLY    EQUIVALENT 

C  BLCCKS    BY    RECORDING    BLOCK    COORDINATE    RANGES    IN     A    MATRIX 

C         3LCCK    BY    THE    USE    OF    RECTANGULAR    COORDINATES 
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C    AS  CUTTING  FUNCTIONS.   THE  CUTTING  COORCINATE  USED  AT 

C    EACH  STEP  IS  CONTAINED  IN  A  VECTOR  IRAD. 

C 

c 

SUBROUTINE    BLCCKR ( VECT ,N ,NN ,M ,MM , IRAD t BLOCK, WKVEC ) 
DIMENSION    VECT(N,M)t BLOCK (NN,MM) ,IRAO(Ni ,UKVEC{M) 
00    10    I=1,MM,2 

BLOCKd,  I)=-1000. 
BLCCK(1,I>1}=1000. 
10  CONTINUE 

DO    100    J=1,N 

00    110    1=1, M 
110  WKVEC(I)=V£CT( J, I) 

00    120    1=1 ,^N 

DO    13C    11=1, M 

IF(WKVEC( II ).LT.BLCCK(  I  ,2*1 I-l)  )G0   TO    120 
IF(i^KVEC(  II). GT.8L0CK(  1,2*11  })G0    TO    120 
130  CGNTINLE 

IBLGCK=I 
GO    TO    150 
120  CONTINUE 

150  JJ=IRAO(J) 

BLIM=WKVEC( JJ) 
DO    160    1  =  1, MM 

BLCCK(J+1,I)=BLCCK( IBLOCK, n 
160  CONTINUE 

BLOCK ( IBL0CK,2*JJ)=6LIM 
3L0CK{ J+1,2*JJ-1)=BLIM 
100         CONTINUE 
RETURN 
END 

C 

C  SUBROUTINE  FOUTR 

C 

C  "  PURPOSE: 

C  TO  CCMPUTE  THE  FOUTZ  STATISTIC  FROM  THE 

C  BLOCKS  DETERMINED  BY  SUBROUTINE  BLOCKR 

C  METHCC  USES  IMSL  ROUTINE  MONCR  TO  EVALUATE 

C  NORMAL  PRCBA6ILITIES  TO  EVALUATE  THE 

C  CLOSED  FORM  EXPRESSION  GIVEN  FOR  D.  THE  FN 

C  STATISTIC  IS  GENERATED  BY  FCUTZ'S  CLOSED 

C  CGMPLTATICNAL  FORMULA. 

C 

C 


SUBROUTINE    FOUTR ( BLOCK ,NN, MM , FN ) 

DIMENSION    BLOCK(NN,MM) ,P(51 ) 

DO    100    1  =  1, NN 

P(  I)  =  1.0 

00    200    J=1,.MM,2 

CALL    MDNOP(BLCCK( I ,J) ,P1J 

CALL    MDNOP{aLCCK( I,J+1 },P2} 

P3=A8S(P2-P1) 

P(  I)  =  P(  I  )*P3 

IF(J. NE. MM-1) GC  TO  200 
200    CONTINUE 
100    CONTINUE 

FN=0.0 

CO  30  0  1=1, NN 

AMAX=1.0/NN-P(I) 

IF(AMAX.LE.O.C)GO  TO  300 

FN=FN+AMAX 
300    CONTINUE 

RETURN 

END 
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